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CN ■ Abstract 



We use characteristic functions to construct a(A;)-multistable measures and integrals, 

Pm . where the measures behave locally like a-stable measures, but with the stability index (x{x) 

varying with time x. This enables us to construct a(x)-multistable processes on IR, that 

»^ ■ is processes whose scaling limit at time x is an a(x)-stable process. We present several 

a . examples of such multistable processes and examine their localisability. 



1 Introduction 
> 

^ I There are several ways of constructing a-stable processes, that is stochastic processes such that 

0^ ■ the finite dimensional distributions of the process at any finite set of m times is an m-dimensional 

a-stable vector, see [lOJ for an full discussion. In this paper we construct a(x) -multistable 
processes, that is processes which look locally like a(jc) -stable processes close to time x in the 

O ! sense that the local scaling limits are a(jc) -stable processes, but where the stability index <x{x) 

varies with time x. 

A number of constructions for multistable processes have been given recently, generalising 

S^ . the constructions of stable processes. One approach is based on Poisson point process [41, 

and another is based on sums of random series [7]. Here we use characteristic functions to 
construct multistable integrals and measures. We show that these multistable measures are 
locally like a-stable measures and may be approximated by sums of many independent a-stable 
measures defined on short intervals with differing a. We then use multistable integrals to define 
multistable processes and give sufficient conditions for processes to be localisable or strongly 
localisable, that is to have a local scaling limit. We give a range of examples of such multistable 
processes. 

2 Definition of a (x) -multistable measure and integral 

Throughout this paper, for given <a <b K^o, the function a : M ^ [a, Z?] will be a Lebesgue 
measurable function that will play the role of a varying stability index. We will work with 
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various linear spaces of measurable functions on M. For < p < oo let 

^p = {/ : / is measurable with ||/||p < 00} where ||/||p =IJ \f{x)\Pdxj ; 

thus II ■ \\p is a quasinorm (i.e. there is a weak triangle inequality ||/ + ^||p < ^(||/||p + ||,?||p) 
for some k> 0) which becomes a norm if 1 < jc < 0°. It is convenient to write 

|/Wr'* = max{|/WIM/W|^}, 

and to define the space of the functions 

^a,b = 7a^7b = {f '■ f'^'^ measurable with / \f{x) I'^'^dx < 00}. 

We also define variable exponent Lebesgue spaces (special cases of Orlicz spaces, see for ex- 
ample [2J) by 



i?^a = {/ : / is measurable with ||/||a < °°} where ||/||a = |x > : f ^ ""^ Jx = l|. 

Then || • ||oc is a quasinorm that reduces to || ■ \\p if <x{x) = ph constant, and is a norm if 1 < a < 
a(x) <b for all x. 

Note that with a < <x{x) <bwe have Ja ^ ^a,b with 

||/||a<Q,,max{||/(x)||„,||/(x)||,}, 

where Ca.t depends only on a and b. 

We define the multistable stochastic integral /(/) of a function / G Ja by specifying the 
finite-dimensional distributions of / as a stochastic process on the space of functions ^a.h and 
then using the Kolmogorov Extension Theorem to show that the process is well-defined. 

Given /i,/2,---,/d G ^a, the following proposition shows that we can define a probabil- 
ity distribution on the vector (/(/i),/(/2), ...,/(/rf)) G W^ by the characteristic function (|)/i,.../-j 
given by (12.11) . The essential point is that a{x) may vary with x. 

Lemma 2.1. Let J G N and /i , /2, ■■■, /d G !Fa, where < a < a{x) < 2 for allx eW. Then 

^fu-fA^U-M = E(exp/f 0,/(/,)) 

j=i 

= exp{-/lE0;//Wr^"^^4 (2-1) 

J j^^ 

for (9i, 02, . . . , Bd) G M , is the characteristic function of a probability distribution on the ran- 
dom vector {I (fi) J {f2),... J (fd))- 



Proof. First, assume that <x{x) is given by the simple function 

m 

a(x) = £aaAiW, (2.2) 

k=\ 

where < a^t < 2 and A^ are disjoint Lebesgue measurable sets with U^^^A^t = ^■ 
For(0i,...,0rf)eR^ 

^^A-l I i e,/;(x) |"« Jxj = exp { - £ / 1 £ e,/, WIa, W r^ j4 

= n^^PI - / 1 Ee,/,wiA,W|"^Jx[. (2.3) 

Now, exp { — / I L/=i 9 j/y (^) Iaj. l-^^) | ''dx] is the characteristic function of the a^t-stable ran- 
dom vector (/(/iIaj), . . . ,/(/rflAj.)), see [[TOll . Hence (12. 3h is the product of the characteristic 
functions of m a^-stable random vectors and so is the characteristic function of a J-dimensional 
random vector given by the independent sum of tty^- stable random vectors. Hence (12.11) is a valid 
characteristic function of a random vector (/(/i), . . .,/(/„)) in the case when <x{x) is a simple 
function (12.21) . 

Now let < a < a(x) < 2 be measurable. Given /i , . . . /rf G Ja write A = {x: Yfj= 1 1 // (x) \ < 
1}. Take a sequence of simple functions {ap(x)}p^ J withO<ap(jc) < 2 suchthatap(;c) -^a{x) 
pointwise almost everywhere; we may assume that (Xp{x) > a{x) if x G A and ap{x) < a{x) if 
X ^ A, for all x and p. Then 

\t^.jM^)r^^^ < maxmax{|0,fM0/}(£|/,W|)"''^^^ 
< maxmax{|0;|Me/}(£|/;W|)"("\ 

J 7=1 

an expression that is integrable since /i, • • ./^ G Ja. By the dominated convergence theorem, 

/ 1 E Qjfj{x)\''"^''^dx ^ / 1 £ e,/,(^) |"(") Jx, (2.4) 

■^ 7=1 -^ 7=1 

and so 

^^p{ [\Y,Qjfj{x)\''''^''~'dx]^cxp{ /l £e,/;W|"(-^)jx|, (2.5) 

■^7=1 -^ 7=1 

as ;? -)■ oo, for all 0i , . . . , 0^ G M. 
For/i,...,/dGiFa, 

/ I I e,/;(x)p(")jx < maxmax{|0/, |0/} £ |/,(^)rW ^ 

-^7=1 ^ 7=1 



as maxy{ |0/| } -^ 0. Thus (O) is continuous at 0. Moreover from (D exp { - / | Z%i Qjfj{x)\ dx} 
is a valid characteristic function of a J-dimensional random vector for all p. Applying Levy's 
continuity theorem to (|2.5I) . there is a probability distribution on the random vector (/(/i ) , /(/2) , . . . , I{fd) ) , 
with characteristic function given by (12.11) . D 

As with a-stable integrals, see [fTOl|, Kolmogorov's extension theorem allows us to define 
a(x)-stable integrals consistently on Ja. 

Theorem 2.2. LetO <a< a{x) < 2. There exists a stochastic process {I{f),f G Ja} with finite- 
dimensional distributions given by A2.1\) . that is with ^/(/i),...,/(/^) = '^fi,...,fdf^^ '^///i, .. .,/rf G 

Proof. For /i, . . . ,/d G ^a it follows from (12.11) that, for any permutation (7l(l),7r(2), ...,7l((i)) 
of (1,2,..., J), we have 

and also that, for any n <d, 

^/„...,/„(ei,...,e„) = (^/„...,/„,...,/,(ei,...,e„,o,...,o). 

Thus the probability distributions given by (12.11) satisfy the consistency conditions for Kol- 
mogorov's Extension Theorem, so, applying this theorem to the space of functions Ja^ there 
is a stochastic process on jFot which we denote by {/(/),/ G jfa}, whose finite-dimensional 
distributions are given by the characteristic functions (12.11) . D 

We call /(/) the a{x)-multistable integral of /. By applying (12.11) with functions (ai/i + 
02/2) J fi, fi and variables 9, —a\Q, — a20 it follows that the multistable integral is linear, that is 
if /i 5 /2 G 'Fa and a\,a2EM., then 

/(ai/i4-a2/2) =ai/(/i)+<32/(/2) a.s. (2.6) 

Let L be Lebesgue measure on M, let E be the Lebesgue measurable subsets of M and 
let "Eq = {A E T, : L(A) < 00} be the sets of finite Lebesgue measure. Let a : M — t- [a,b] be 
measurable where < a < b <2. Analogously to [fTOl Section 3.3] for a-stable measures, we 
define the a{x) -multistable random measure M by 

M{A)=I{1a) (2.7) 

for A G T.Q, where 1a is the indicator function of the set A; thus M{A) is a random variable for 
each A G Eq- 

It is natural to write 

lf{x)dM{x):=I{f), /G5a, (2.8) 

since there are many analogues to usual integration with respect to a measure. Clearly, linearity 
of this integral is a restatement of (12.61) . and 

lAix)dM{x)=M{A). 



With this notation the characteristic function (12.11 ) may be written 

E(exp/{£0,|/,(x)JM(x)})=exp|-||£0,/,W|"^")jx} (2.9) 



for fj E 7a- For the random measures, taking fj = 1a with Ay G Eo. 



E(exp/{ £ 0,M(A,) }) = exp I - / I £ G^U,.! 



x)| (i;c 



(2.10) 



We may estimate the moments of an a(jc)-multistable integral in terms of the norm || ■ ll^. 

Proposition 2.3. Let < a < (l{x) <b<2 and let g G Ja- Then there is a number ci depending 
only on a and b such that for allX> Q 



P 



I g{x)dMa{x)\>'k\ <c,j 



fW 



X 



a{x) 



dx. 



(2.11) 



Moreover, ifO < p < inf;(^£Ka(x) there is a number C2 depending only on p and b such that 



E 



g{x)dMa{x) 



<C2||^||a- 



(2.12) 



Proof. A simple calculation using distribution functions (see [ill p.47]) gives 



P 



g{x)dMa{x] 



J 


< 


2 J-2/X 




= 


2 J-2/X 




< 


X /■2/^ 
2 J-2/X 



1 - E^exp {iQjg{x)dMa(^,))UdQ 
('l-exp(- /"|0g(jc)rWjjc))j0 



|0|«W|^(x)|''Wjx)rfe 



< ci 



I 



^w 



a{x) 



dx. 



Assuming as we may that ci > 1 and writing Xq = \\g\\a > for the number such that 



/Xo"^'^|^(jc)|«Wjjc = 1, we have 



E 



g{x)dMa{x] 



p I XP-^F 

'0 



>X]dX 



g{x)dMaix) 

< cip Tv^-Vinjl, /"x-'^WlgWrWjjcjjX 

< cip xp-^dx+cip ;i^-i-''('^)|g(jc)|"Wjxjx 

Jo J JXo 

< c,Xl + c,Xl j X-Q''^'^\g{xT^'Ux 



= ClUWa- 



u 



Recall that a random measure M on M is independent scattered if M(Ai),M(A2), ...,M(Arf) 
are independent whenever A i^2v:^/t ^ "^Q ^e pairwise disjoint, and is o-additive if whenever 
Ai,A2,...G £0 are disjoint and [JJ^jAy G £0 then almost surely 

7=1 )=i 

taking an independent sum. 

Theorem 2.4. The a{x)-multistable measure Ma is independent scattered and o-additive. 

Proof. This is a slight variant of |[TOl Section 3.3]. Let Ai,A2,.--^A: ^ "^o be pairwise disjoint. 
Then using (ITTOl) 



.)})■ 



E(exp{/£0;Ma(A;)}) =nexp|- /|0,lA,Wr^"^j4 = n^(e'^P{''9i^«(^ 

so Ma(Ai),Ma(A2), ...,Ma(Arf) are independent, and Ma is independent scattered. 
If Ai,A2,...,Aa: G £0 is a finite collection of disjoint sets, using (12.71) and (12.61) . 

k k k k 

Ma{[JAj) =I{Ua,) =^(I1a,) = I:/(1a,) = l^MaiAj). 

7=1 7=1 7=1 7=1 

For a countable family of disjoint sets Ai,A2, ... G £0 with 5 = 1J7=1^7 ^ "^O: so that 5 = 
\j)^iAj U [ U7=fc+i^7) ' it follows from above that 

k ao k '=° 

Ma{B)=Ma(\jAj)+Ma( [J Aj) = Y^Ma{Aj) +Ma( [J Aj). (2.13) 

7=1 j=k+l 7=1 ^j=k+l 

Since lim^^oo L ((J7=A:+i^7) = ^ and a{x) G [a, Z?], for each G IR 



limEfexp/{0Ma( M Aj)}) = limexpj - / |01 



j=k+l 



L)j^,+,Ajr'--' ^-1. 



d 



SO Ma{\J7=i^^iAj) — i- as fc — !■ 00 by Levy's Continuity Theorem. 

By (1233]) we get Ma{B) - Z^iMaiAj) A and so Ma{B) - E5=iMa(A;) A as A: ^ 

00. Thus \imk^'^'L)=i Ma{Aj) ~ ^a{B), and, since the summands Ma{Aj) are independent, 
this implies convergence almost surely, by a theorem of Kolmogorov, see f6l. Thus Ma is 
a- additive. D 

Next we obtain conditions for convergence of a sequence of multistable measures with dif- 
ferent multistable indexes. 



Proposition 2.5. Let a„(x), a(x) be Lebesgue measurable with < a < a„(x),a(x) <b<2 
for all X &W. Let Ma„,Ma be the associated (X„{x)-multistable and a{x)-multistable measures 

characterised by f |2.i0|) . Suppose a„ {x) — )■ a{x) for almost all xeM.. Then Ma„ — > Ma asn-^oo, 
that is for all m G N andAi,A2, ■■■,Af„ G £0- 

(Ma,.(Ai),Mcx„(A2),...,Ma„(A„)) A (Mcx(Ai),Ma(A2),...,Ma(A™)). 
Proof Let Ai, A2,."4ra G Eo- Then for all n and all x G M 

m 

7=1 

where A = UyLi^j e Eo and c = max { (^Li |e;|)^ (E"Li |0;i)*}- Since flA{x)dx < 00, the 
dominated convergence theorem implies that 

Jimexp(-|| f;0;U/x)|""^"^J^) =exp(-|| £0,1a,.(x)|"'("^Jx), 



7=1 

SO from (12.101) . 



7=1 



E(exp/{ £ BjMaMj)}) ^ E(exp/{ £ e,Ma(A,)}) 

7=1 7=1 

as n — > 00. By Levy's continuity theorem Ma„ — )■ Ma. D 

To get a feel for a(jc)-mutistable measures, we show that, for a continuous a{x), the a{x)- 
multistable measure M may be approximated by random measures that are the sum of many 
independent a-stable measures defined on short intervals. 

Assume that a : M — )■ [a, Z?] C (0, 2] is continuous and let Ma be the a(x)-multistable measure 
on the sets "Eq. We now use the same procedure but with piecewise constant functions a„(jc) to 
obtain approximating measures Ma„- 

For each n let a„ : M — )■ [a, b] C (0, 2) be given by 

Unix) = a{r2-") ifxe [r2-'\ (r+ 1)2-") for r G Z 

and let Ma„ be the a„-multistable measure obtained from an{x) as above, so in particular Ma„ 
has finite-dimensional distributions given by the characteristic function 

E(exp/{ f QjMa„{Aj)}) = exp { - / | £ 0,1^. Wf-^^^Jx). 
7=1 ■^ 7=1 

It follows from Theorem [2i4l that each Ma„ is independent scattered and o-additive. 



Theorem 2.6. Let 0<a<b<2 and a : M — )■ [a, Z?] be continuous. LetM„,- denote the restriction 
of CL{r2^") -stable measure to the interval [r2^", (r + 1)2^")), that is 

M„,,(A)=M„(,2-„)(An[r2-",(r+l)2-"))=M„„(,)(An[r2-«,(r+l)2-«)), 

where M^ty2->^) '■^ (X,{r2^") -stable measure. Then Ma„ is a random measure given by the inde- 
pendent sum of random measures 

reZ 

almost surely for A E "Eq. Moreover Ma„ — )■ Ma as n ^ oo. 

Proof. Since M^^ is independent scattered, we have that for each A E 'E 

Ma„(An[r2-",(r+l)2-'')) =M„,,(A) 

are independent for distinct r. 

Let A E Eo- Since Ma„ is o-additive, 

M„(A) = Ma„iA) 

= Ma„{[JAn[r2-\{r+l)2-")) 

= £Mo.„(An[r2-",(r+l)2-")) 

reZ 

where the summands are independent. 

For each n we have a„(;c) = a(r2^") for all x E [r2^'\ {r+ 1)2^"). Since a{x) is assumed 

fdd 



continuous, we have lim„^oo cx„ (x) = a{x) for all x. Thus by Theorem 12.51 Ma„ — )■ Ma as 

n — 7- oo. n 

One would expect an a(x)-multistable measure to 'look like' an a(a)-stable measure in a 
small interval around u. We now make this idea precise. 

For M G M, r > 0, let Tu^r : M — t- M be the scaling map. Tutrix) = {x — u)/r. This induces a 
mapping T*^. on random integrals and measures, given by 

f f{x)d[TlMa){x)= f f[^)dMa{x)^l[f[-^)). 
In particular 

(r*,Ma)(A) =M„(r,-i(A)) =/(i^-,(^)) 

for AG En bv (l2Jl). 

We show that scaling an a(x)-multistable random measure about a point u yields the <l{u)- 
stable measure Ma[u) ■ 



Theorem 2.7. Let a : M — )■ [a,Z?] C (0,2] be continuous with 

|a(x + r)-a(jc)| =o(l/logr) (2.14) 

uniformly on bounded intervals and let m G M. Then for all functions /i , . . • , /j G Ja,b with 
compact support, the vectors 



4 (|/iWJM„(,)W,...,|/dWJM„(,)W) (2.15) 
a* r — )■ 0. In particular, 

^-l/a{u) ((7;#^Ma)(Ai), . . . , (r*,Ma)(Arf)) A (M„(„)(Ai), . . . Ma(u){Ad)) (2.16) 

a* r — > 0, /or all bounded sets Ai , . . . , A^ G Eq- 

Proo/ Let /i,/2, ...,/m G J^^^ be functions with compact support, say in [— zo^zo]- Let 0^ G M, 
7 = 1, 2, ..., m, and consider the characteristic functions. 

E 



(exp/£0,r-i/'^W//,WJ(r*^a)O 

^ 7=1 •' 

= E(exp/|;0,r-i/«(")|/,(^)jM„( 
= exp(-/|f0,.-VaH^^.(^)|«W,,) 

= exp(-/|£0,r-iM")/,(z)r(^^+"Vjz) 

= exp(- /|£0,/;(z)p(''^+"Vi-'^('-^+")M")jz), (2.17) 

on writing (.jc — u)/r = z. 

From condition (|2.14l) it is easy to see that Um^^o''^ " ^^^" = 1 uniformly for z G 

[— zo,zo]j and also lim^^oOc(''2 + ") = Oc(") uniformly for all z G [— zo^zo] since a is continuous. 
Noting that there is a constant c such that for r sufficiently small, 

m in 

I E 0,/;(z)p('"^+"Vi-«('-^+")/«M < c £ \fj{z)\"'\ 

7=1 7=1 

for z G [— zo:Zo] and fj G ^a,/?. the dominated convergence theorem gives 

limexp(- / I £0,/;(z)p(^^+"Vi-'^('-^+")/«Mjz) =exp(- / | £ 0;/;W|"(")jx), 



soby dTTTl) 

limE(exp/£0,r-iM")y /;WJ(r*^cx)W) =E(exp/£0;y /,WJM„(,)W). 
Levy's continuity theorem now implies (12.151) and (I2.16I) . D 

3 Multistable processes and localisability 

In this section we introduce processes defined by multistable integrals, and in particular consider 
their local form, with the aim of constructing processes with a prescribed local form. Thus, 
given a : M — )• [a, Z?] C (0, 2], we write 

Y{t) = j f{t,x)dMa{x), (3.1) 

for ? G M and / G Ja,/„ where the integrals are with respect to an a(x) -multistable measure M^ 
as in (|2.8I) . By (12.91 ), for each {t\,t2, ■■■,td) G W^ , the characteristic function of the random 

vector (y(?i),y(?2),...,y(^d)) is 

E(exp/£0;y(?;)) = E{expiY, f f{tj,x)dMa{x)) 
= exp (^ - y I £ Qjf{tj,x) I ''^'''^dx 



forall(0i,02,...,erf)GM^'. 

First we give conditions for Y to have a continuous version. 

Proposition 3.1. Let a : M — t- [a,Z?] C (1,2] be measurable and suppose fit,-) G la for all 
t E^. Let Y be given by ( I3.il) . Suppose that there exists \/a <V[ < \, such that for each 
bounded interval I we can find c > such that 

||/(f, ■)-/(", ■)||a<c|f-wn {t,u(^I). (3.2) 

Then Y has a continuous version satisfying an a.s. ^-Holder condition on each bounded interval 
for a// < P < (r\a — \)/a. In particular, ( 13.21) holds if 

f\f{t,x)-f{u,x)\''^''Ux<ci\t-up (t^uEl), (3.3) 

a form that may be easier to check in practice. 

Proof. Take p such that l/r\< p < a.By Proposition [23] 

E{\Y{t)-Y{u)\P)=E(\l{f{t,x)-f{u,x))dMa{x)\')<C2\\fit,-)-fiur)fa<c^^^^ 

The conclusion follows from Kolmogorov's continuity theorem, see [9, Theorem 25.2]. D 
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Recall that a stochastic process Y is localisable at a point if it has a unique non-trivial scaling 
limit, formally Y = {Y{t) : ? G M} is h-localisable at u with local form or tangent process 

Y:, = {Y:,it):teR}i{ 



Y{u + rt)-Y{u) m^, 



(3.4) 



as r — 7- 0. If y and 7^' have versions in C(IR) (the space of continuous functions on IR) and 
convergence in (|3.4I) occurs in distribution with respect to the metric of uniform convergence 
on bounded intervals we say that Y is strongly localisable. For the simplest example, a self- 
similar process with stationary increments Y is localisable at all u with local form Y', = Y and is 
strongly localisable if it has a version in C(IR) . In general there are considerable restrictions on 
the possible local forms, see S. 

We call a stochastic process {Y{t), ? G M} multistable if for almost all u, Y is localisable at 
u with Fj' an a-stable process for some a = <x{u), where < <x{u) < 2. Various constructions 
of multistable processes are given in ||4l[51l3. 

For a stochastic process Y, it is natural to ask under what conditions Y is localisable. The 
following theorem, which is a multistable analogue of |l5l Proposition 2.1], gives a sufficient 
condition. 



Theorem 3.2. Let 



Y(t) 



I f{t,x)dMa{ 



(3.5) 



where Ma is an a,{x) -multistable measure for continuous a : M — )■ [a,Z7] C (0,2]. Assume that 
fit, .) G Ja,bfor all t and 



lim / 



f{u + rt,u + rz) - f{u, u + rz) 



Ji—i/a{u+rz) 



-K^z] 



a,b 



dz = 



(3.6) 



for a jointly measurable function h{t,z) with h{t, .) G fa.bfor all t. Then Y is h-localisable at u 
with local form 

y:=||/z(^z)JM„(,)(z):rGM} (3.7) 

where M^iu) i^ (x{u)-stable measure. 

If in addition, there exists r\> 1/a such that for each bounded interval I we can find c> 

such that 

f{u + rt,-)-f{u + rv,- 



<c\t-v\'^ (t.vel) 



(3.8) 



for all sufficiently small r>0, then Y is strongly localisable at u. Condition 63.8\l is implied by 

f{u + rt,u + rz) - f{u + rv,u + rz) «("+'-z) 



fh~l/a{u+rz) 

which can be more convenient to use in practice. 

To prove Theorem 13. 2[ we need some convergence estimates 



dz<Ci\t-vp (t.vel) 



(3.9) 
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Lemma 3.3. Let < a < b. There is a constant c that depends only on a and b such that, for 
all measurable a : M — )■ [a,Z?] and g,k G 'Ta,b' 

^ (w ; II II 7 i|max{0,a— 1} , n , na , n , n II , ||max{0,fc— 1} , n i \\b\ /o i rv\ 

<c[\\g-k\\a\\k\\a +||^-fc||^+||^-^|U||^||^ ^ + \\g-k\\iy (3.10) 

Proof. If < fl < a{x) < Z? < 1 for all jc G M, then 

On the other hand, if 1 < a < <x{x) < b for all x G M, then by the mean value theorem there 
exists < X{x) < 1 such that 



la— 1 
\b-\ 



b-l 
b 



< b\\g{x)\-\k{x)\\\\k{x)\ + \g{x)-k{x) 
+b\\g{x)\-\k{x)\\\\k{x)\ + \g{x)-k{x) 
Integrating and using Holder's inequality gives 

which gives (13.101) in this case. 

In general, for < a < <x{x) < b, letting A = {x: a < a{x) < 1}, inequality (13.101 ) holds for 
glA and klA and also for gl^\A and fclR\A' ^^'^ combining these cases we get (13.101 ) for g and k 
for an appropriate c. D 

We require the following Corollary. 13. 2[ 

Corollary 3.4. LetO<a<b and g:R+ xR^R* with g{r, .) G Ja,bfor all r>0. Letke fa,b 
and let^:R^ [a, Z?] be continuous at 0. If 



\im f\g{r,z)-k{z)r'''dz = 0, 



then 



limf \g{r,z)f'-'Uz = J |/r(z)|P(OW 
Proof By dlH]) and Lemma [33] 

lim f\g{r,z)\^^'-'Uz- /"|fc(z)|P('-^W 
Since k G !Fa.b, the dominated convergence theorem gives 

lim f\k{z)f'^Uz- f\k{z)\^^^Uz 

r^O J J 

and (|3.12l) follows on combining these two limits. 



(3.11) 
(3.12) 



0. 



0, 



D 
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We can now complete the proof of Theorem |3.2[ 



Proof of Theorem 13.21 Fix m G M. We consider the characteristic function of the finite- 
dimensional distributions of r~^{Y{u-\-rt) —Y{u)). Let 0^ G M and tj G M for j = l,2,...,m. 
Then, using ([33]) and (IZ9l) . 

171 

^Ux^iY^Qjr-^{Y{u + rtj)~Y{u))\ (3.13) 

= E(exp/£0^r"'' f {f{u + rtj.x) - f{u,x))dM{x)) 
= oxp\ - [ \f^Qjr-''{f{u + rtj,x)-f{u,x))\''^'Ux] 



(3.14) 



after setting x = rz + m. 
Defining 

Z{t)= Jh{t,z)dMa(u)iz), 
its finite-dimensional distributions are given by the characteristic function 



E(exp/£0jZ(f^)) =expj - / \Y,Qjh{tj,z)p"'dz\. (3.15) 

We now use Corollary 13. 4[ taking 

. _y. f{u + rtj,rz + u)-f{u,rz + u) 
7=1 r ' y 

m 

kiz) = l^Qjh{tj,z), 

7=1 

and 

(3(jc) = a(M+.x:). 

Then 

||g(r,z)-fc(z)rJz^O, 

as r — 7- 0, using (|3.6l) and the quasi norm properties of || • ||a and || ■ ||^. Thus by Corollary 13.41 

/m I \ r "^ t \ 

I £0^.,-/^+lM.^+")(/(^ + rt„rz + M)-/(i<,rz + M))p(''^+")jz^ / I £0,;Mo-,z)r^"^rfz, 
7=1 •' 7=1 
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as r — )■ 0. 

Since the exponential function is continuous, (13.141) . and hence (13.131) . is convergent to 

(I3l5l) as r ^ for all (Bi, . . . ,0„). By Levy's Continuity Theorem, r-'\Y{u + rt) - Y{u)) ^ 
Z{t) as r — )■ 0, noting that (13.151) is a characteristic function. Thus Y is /z-localisable with local 
form y^' given by (13.71) . 

Finally, if (13.81) holds then by Proposition [231 for < p < a. 



< C2 



f{u + rt,-)-f{u-rv,-) P 

a 



yh 



< cslr-v^p. 



By choosing p such that l/r\< p <a, Kolmogorov's continuity theorem, see DH Theorem 25.2], 
implies that, for each < P < {r[a — l)/a and each bounded interval /, the process 7^ satisfies 
an a.s. Holder condition 

where the random constants behave uniformly in r, i.e, supQ^,.<;,.yP(C^ > m) — t- as m — )■ oo. 
Thus for all £, x > there exists 5 > such that 

limsupP( sup \Yr{t) - Y,.{v)\ > xj < £. 

In other words, the Yy are strongly stochastically equicontinuous on / which, along with conver- 
gence of the finite dimensional distributions, implies that Yr converges to Y' in distribution on 
the space of continuous functions with the metric of convergence on bounded intervals, see |[I1 
Theorem 8.2] or[8, Theorem 10.2] . □ 

4 Examples 

We give a number of examples to illustrate Theorem 13.21 Some of these are considered in 
[|4l[5l|71 using alternative definitions of multistable processes. 
It is convenient to make the convention that 

if V < M in the following examples. 

Example 4.1. Weighted multistable Levy motion. 
Let 



Y{t) = jw{x)lioj]{x)dMa{^){x), 
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where CL:M.-^ [a,2\ is continuous and a> 0, and w :M.—>'Ris continuous. Let u EMbe such 
that w{u) 7^ and suppose that asv^u, 

\a{u) -a(v)| =o(l/|log|M-v||). (4.1) 

Then Y is \ / (l{u)-localisable at u with local form 

K = { Jw{u)lio^t]{z)dMa(u){z), r e m} = w{u)La(u), 

where Lot(-„) is a a{u)-stable Levy motion 

Proof. Take f{t,x) = w(x)l[o,/](x) and h{t,z) = w(m)1[o,/](z). Condition (14.11) ensures that 
j,i/a{u)-i/a{u+rz) _v. 2 as r — )■ uniformly for z G [0,f] which is needed to ensure that (|3.6I) 
holds. Then Theorem 13.21 gives the conclusion. D 

Next we consider multistable reverse Ornstein-Uhlenbeck motion. Notice that in the multi- 
stable case, we get a curious restriction on the range of a. 

Example 4.2. Multistable reverse Ornstein-Uhlenbeck motion. 
Let 

exp(-X(x-0)^M„(,)(x), (4.2) 

where a: M — )• [a,b] C (1,2] is continuous with 1 < \/b < a < b <2. Let u eW and suppose 
that asv-^u, 

\a{u) -a(v)| =o(l/|log|M-v||). (4.3) 

Then Y is l/(l{u)-localisable at u with local form 

y:={|-l(o,,)(z)JM„(,)(z), teR]. (4.4) 

Proof. We take f{t,x) = exp{—'k{x — t))l[f^^^{x) and h{t^z) = — l[o.f)(z) in Theorem [X2l After 
a little simplification, 

f{u + rt,u + rz) -f{u,u + rz) .. , "'^ 



/■ f[u + rt,u + rz)-f[u,u + rz) 
J f.l/a{u)-l/a{u+rz) ^^'^ 

r\t\ -exp{-'krz)l[oj)iz) a,b /■- exp(-Xrz)(exp(Xr?) - l)\yt,^){z 

~ J-\t\ r^/a{u)-i/a{u+rz) + Ho,tM) ^^ J^^^ f.l/a{u)-l/a{u+rz) 

The first integral converges to 0, noting that j-i/«(")-i/«("+''z) _!. i as r — ;■ 0, uniformly on z G 
[— f^,;^]. The second integral is bounded by 



dz 

a,b 



dz. 



/'OO 
I r- 1/^^+1/* exp(-Xrz)(exp(XrO - \)\'dz 
A 

/'OO 
|exp(-?irz)(exp(Xrt)-l)r'^^z 
A 

poo 

< ciri-^/"|exp(XrO-ir' / \exp{-Xarz)\dz 

J \t\ 



< C2r^-''/''{Xr\t\yexp{-Xar\t\){Xray^ 
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for fixed t, where c\, C2 and C3 are independent of r < 1 . Since a — b/a > the second integral 
converges to 0, so the conclusion follows from Theorem 13 .21 D 

The next example is linear fractional multistable motion. Recall from IfTOl that asymmetric 
linear fractional a- stable motion is given by 

/oo 
Pa,/.(^+,^",^-^)^MaW (4.5) 

-00 

where t,b~^,b~ eM., and 

andM(x is a-stable random measure (0 < a < 2). By convention, if /z — 1/a = 0, we take 

if ? > 0, and 

pa,h{b^ :b^ ,t,x) = -{b^ - b^)lltfi]ix) 

ift <0. Then (14.51) is an a-stable process. 

For a multistable version let a : M — )• [a^b] C (0,2) be continuous. We define linear frac- 
tional a{x) -multistable motion by 

/oo 
Pa{x),h {b^ , b^ , t,x)dMa^^^ (x) (4.6) 

-co 

where t eR,b^,b^ eM, and 

p,^^^^,ib\b-J^x) = b^it-x)'-'/^^^^-{-x)i-'/''^'^) 

where Ma{x) is a(x) -multistable random measure. 

It may be checked directly that if ? G Mand l/a—l/b<h< \ + l/b— l/a\he,n^^L\j^{b^,b^,t, 
fa,b SO that (14.61) is well-defined. 

We show that linear fractional multistable motion has linear stable motion as its local form. 
We consider the case when b^ = 1 and b^ = 0, the argument is similar for other b^ and b^ . 

Proposition 4.3. Linear fractional multistable motion. 
Let 

Y{t) = /(r-;c)l-i/"«-(-^)Ti/«WjM,w(;c) 

= y Pa(x),/,(l,0,^^)^^a(x)(-«) 
= -^aW, A, 1,0(0' 
where a: M — t- [a,b] C (0,2) is continuous. If 

\/a-\/b<h<l + l/b-\/a, (4.7) 
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then Y is h-localisable at each w G M with local form 

= ■^a(«),/!,l,o(0- 

Furthermore, if \/a < h < 1 + 1/b — 1/a, then Y has a continuous version and is strongly 
localisable at each u &M.. 

Proof We take f{t,x) = (f -x)'!"^/"^""^ - (-jc)'!"^/''^'^) e 5a,&, given (gj]), and h{t,z) = [t - 
^^h-i/a{u) _ ^_^^h-i/a.(u) .^ xheoremm Then 



/ 



f{u + rt, M + rz) - f{u, u + rz) 



r 



.h-i/a{u+rz) 



-Khz] 



a,b 



/ \h-l/a{u+rz) _ f _ \h-l/a{u+rz) _ j \h-l/a{u) , _ ■.h-l/a{u) "' 

I' ^J+ ~\~Z)+ ~1'~ZJ_|_ +[—z)^ 



dz 

dz. 



This integral converges to as r — )■ 0. This may be established by breaking the range of inte- 
gration in the parts: |z| < 5, |z — f| < 5, |z| > M and A = {z : 5 < |z| < M and 5 < |z — f|}. By 
choosing sufficiently small 5 and large M, the integral over the first three parts can be made arbi- 
trarily small, uniformly as r — )• 0. The integrand converges to pointwise on A and the bounded 
convergence theorem gives the integral over A convergent to 0. The conclusion follows from 
Theorem l3.2[ 

Finally, if 1/a < /z < 1 it is easily checked by a routine integral estimate (lit > u splitting 
the resulting integral at u) that (|3.3I) holds so Y has a continuous version, and similarly that (|3.9I) 
if 1/a < T] < h, so Y is strongly localisable by Theorem 13 .21 D 
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